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Abstract. The discovery of millisecond pulsations from neutron stars in low mass X-ray binary 
(LMXB) systems has substantiated the theoretical prediction that links millisecond radio pul- 
sars (MSRPs) and LMXBs. Since then, the process that produces millisecond radio pulsars from 
LMXBs, followed by spin-down due to dipole radiation has been conceived as the 'standard evo- 
lution' of millisecond pulsars. However, the question whether all the observed millisecond radio 
pulsars could be produced by LMXBs has not been quantitatively addressed until now. 

The standard evolutionary process produces millisecond pulsars with periods (P) and spin-down 
rates (P) that are not entirely independent. The possible P — P values that millisecond radio pulsars 
can attain are jointly constrained. In order to test whether the observed millisecond radio pulsars are 
the unequivocal descendants of millisecond X-ray pulsars (MSXP), we have produced a probability 
map that represents the expected distribution of millisecond radio pulsars for the standard model. 
We show with more than 95% confidence that the fastest spinning millisecond radio pulsars with 
high magnetic fields, e.g. PSR B1937H-21, cannot be produced by the observed millisecond X-ray 
pulsars within the framework of the standard modelQEI 
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INTRODUCTION 

Normal pulsars (P ~ few seconds, B ~ lO^^G) in low mass X-ray binary (LMXB) 
systems are believed to potentially spin-up to reach milliseconds rotational periods by 
acquiring angular momentum from their companion fUH^]- 

The millisecond modulations in the X-ray signal observed from some of these actively 
accreting neutron stars and bursting sources are likely to be signatures indicative of their 
spin periods. There are about ~20 nuclear or accretion powered (see Table 1) millisec- 
ond X-ray pulsars (MSXPs) which are thought to be the progenitors of millisecond radio 



pulsars (MSRPs) [|39ll . These millisecond X-ray pulsars may become observable in ra- 
dio wavelengths once accretion ceases, or the column density of the plasma from the 
fossil disk around the neutron star becomes thin enough to allow vacuum gap formation 
that leads to the production of coherent radio emission. Towards the end of the secular 
LMXB evolution, as accretion rates fall below a critical value above which detection 
presumably may be hampered due to absorption or dispersion [|37|] . the neutron star can 
re- appear as a millisecond radio pulsar. 



1 For details see Kiziltan & Thorsett (2009) 

^ High resolution full color figures available at URL: |http : // w w w.kiziltan . org/research.htrnl 



Although the connection between LMXBs and MSRPs has been significantly 
strengthened after the discovery of quasi-periodic kHz oscillations and X-ray pulsations 
in some transient X-ray sources 31, 12, [13] > no radio pulsations from millisecond 
X-ray pulsars have been detected so f^rljr 

Accretion physics [fill] predicts an equilibrium period (Peq ~ [m/MEdd]^^^'^) of 
the neutron star that is accreting in the LMXB system which is constrained by the 
accretion rate (m). This equilibrium period combined with the dominant mechanism for 
energy loss delineates the subsequent kinematics of the spun-up millisecond pulsar. The 

magnetic dipole model then implies a "spin-up region" (P ~ Pq^^) [seef^] on which the 
recycled neutron stars will be reborn as millisecond radio pulsars. At the end of the active 
phase, millisecond X-ray pulsars accreting with m and spinning with Peq, presumably 
transition into a millisecond radio pulsar with an initial spin period of Pq ~ Pgq. 

The energy loss due to magnetic dipole radiation after the onset or radio emission 
is the main process that drives the evolution for millisecond radio pulsars [see 



2M- Alternative energy loss mechanisms such as multipole radiation or gravitational 



wave emission, especially during the initial phases of the reborn millisecond pulsars, 
have been suggested by several authors [f22L [3] but have yet to be observationally 
corroborated. Advanced Laser Interferometer Gravitational Wave Observatory (LIGO) 
will be able to probe the frequency space at which millisecond pulsars are expected to 
radiate gravitational waves, thereby putting stringent constraints on the micro physics of 
millisecond pulsars. 

The advances in radio observations, increased sky coverage with deep exposures of 
current surveys combined with robust post-bayesian statistical techniques that incorpo- 
rate minimal assumptions, give us unprecedented predictive power on the joint period- 
spindown (P — P) and implied magnetic field (B) distributions. 

Here, we attempt to go beyond phenomenological arguments and test whether mil- 
lisecond X-ray pulsars can produce the characteristics of the observed millisecond radio 
pulsars within the framework of the standard model [i3i, and references therein]. 



CONSTRAINING MILLISECOND PULSAR EVOLUTION 

In order to statistically constrain millisecond pulsar evolution, we need to describe the 
parameters that drive the evolution in quantifiable terms. The evolution of millisecond 
pulsars can be consistently described in terms of i) the equilibrium period distribution 
(D) of millisecond X-ray pulsars at the end of the LMXB evolution ii) the mass accretion 
rates (M) of the progenitor population during the recycling process iii) Galactic birth 
rates (R), and iv) the dominant energy loss mechanism after the onset of radio emission. 

We devise a semi-analytical evolution function S' to parametrize the evolution of 
millisecond pulsars after the accretion phase, which can be described in closed form 
as: 

i=0 1=0 

where i^Qi^ is the probability distribution function. The shape parameters a and /3 



define the distributions (i.e., D,M,R for k=l,2,3) for the Beta function^ [9] inferred 
from observations at each Monte-Carlo realization "r". 

The evolution function S' is built by randomly choosing initiation seeds from a 
period distribution D, which is then convolved via the standard model to consequently 
sample the P — P parameter space. For the observed millisecond X-ray pulsars, the 
period distribution which seeds will be randomly chosen from is the observed Pmsxp 
distribution (Table 1). We uniquely construct a "relaxed multidimensional Kolmogorov- 
Smirnov (K-S) filter" (Fig. 1) to check population consistencies by calculating the 2D 
K-S M] probabilities {Pidk-s) between observed millisecond radio pulsars and the 
synthetic population that is formed by these properly evolved progenitor seeds. The 
filtering is reiterated for each realization to obtain synthetic populations with consistent 
distributions. The filtered (D,M,i?) values are then convolved with standard dipole 
energy loss in order to construct the probability distribution function {S^Si^) which 
is the integrated form of the proper {P,P) pairs. 

Nominally any P2DK-S > 0.2 value would imply consistent populations in a 2D K-S 
test. By allowing 0.005 < P2DK-S < 0-2 with lower fractions (see Fig. 1), we oversample 
outliers to compensate for possible statistical fluctuations and contaminations. A peak 
sampling rate around the nominal acceptance value of P2DK s ~ 0.2 is the most optimal 
scheme that prevents strong biases due to over or under- sampling. The main goal for 
oversampling outliers and relaxing the K-S filter is to test whether the standard model 
can at least marginally produce very fast millisecond pulsars with relatively high mag- 
netic fields like PSR B 1937-1-21. 

The predictive significance of the P — P distribution for the probability map (Fig. 2) is 
obtained from a Monte-Carlo run with r = 10^ valid realizations that produce consistent 
synthetic samples. Whilst sampling the P — P space, no assumptions were made regard- 
ing the progenitor period distribution (D), the accretion (M), or the Galactic birth {R) 
rates. The filter (Fig.l) is implicitly driven by the observed millisecond radio pulsars. 

Figure 2 shows the expected P — P distribution for the standard model assuming 
that millisecond radio pulsars have evolved from a progenitor population similar to 
the observed millisecond X-ray pulsars. We do not include millisecond radio pulsars 
in globular clusters because the P — P values in these cases may not necessarily be 
the sole imprint of the binary evolution, but can be significantly changed by possible 
gravitational interactions due to the crowded field. To explore the extend of the effects 
of an unevenly sampled progenitor population, we also show the region in the P — P 
space that is sensitive to alternative Pmsxp distributions. The probability map is overlaid 
with the observed millisecond radio pulsars. 



^ Beta functions are commonly preferred in Bayesian statistics as the least restrictive and most flexible 
prior distributions. It can take the form of an uninformative (e.g. uniform) prior, a monotonic line, concave, 
convex, unimodal (e.g. normal) or any extreme combinations of these shapes. 



CONCLUSIONS 



• Young millisecond pulsars with higher magnetic fields (e.g. PSR B1937+21) are 
inconsistent with the standard model. The fraction of the observed young/old mil- 
lisecond radio pulsars with high B fields is higher than what the standard model 
predicts by several orders of magnitude. 

• The fastest spinning millisecond pulsars, in particular PSR B 1937+21, may origi- 
nate from a different evolutionary channel. 

• The standard evolutionary model is able to produce the general demographics of 
older millisecond radio pulsars only. 

• Accretion rates that millisecond radio pulsars have experienced during their ac- 
cretion phase deduced from observed P — P values, combined with the observed 
millisecond X-ray pulsar period distribution produces mostly older millisecond ra- 
dio pulsars, including millisecond radio pulsars with spin-down ages Xc > 10^*^ yrs. 
Figure 2 shows clearly that the apparent enigma of millisecond pulsars with spin- 
down ages older than the age of the Galaxy is mainly a manifestation of very low 
accretion rates during the late stages of the LMXB evolution. 

• No physically motivated P distribution has been able to produce the whole mil- 
lisecond radio pulsar population consistently. No millisecond X-ray pulsar period 
distribution could mimic the observed relative ratios of young/old pulsars with high 
B fields. 

• It is necessary to posit the existence of a separate class of progenitors, most likely 
with a different distribution of magnetic fields, accretion rates and equilibrium spin 
periods, presumably among the LMXBs that have not been revealed as millisecond 
X-ray pulsars. 

• A millisecond X-ray pulsar period distribution that has sharp multimodal features 
coupled with non-standard energy loss mechanisms may be able to reconcile for 
the joint P — P distribution of millisecond pulsars. 
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TABLE 1. Accretion and nuclear powered 
pulsars: The millisecond pulsar progenitor seeds 
used to construct the cumulative synthetic mil- 
lisecond radio pulsar population for the observed 
Pmsxp- 



Vspin [Hz] 


Pulsar 


Reference 


619 


4U 1608-52 


[16] 


601 


SAXJ1750.8-2900 


[18] 


598 


IGR J0029 1+5934 


L30] 


589 


X 1743-29 


[35] 


581 


4U 1636-53 


[42] 


567 


X 1658-298 


m 


549 


Aql X-1 


[4r] 


524 


KS 1731-260 


[Ml 


435 


XTEJ1751-305 


[31J 


410 


SAXJ1748.9-2021 


\M 


401 


SAX J1808.4-3658 


[39,7] 


377 


HETE J 1900. 1-2455 


[20] 


363 


4U 1728-34 


{36} 


330 


4U 1702-429 


\2E 


314 


XTEJ1814-338 


[2^ 


270 


4U 1916-05 


[la 


191 


XTEJ1807 .4-294 


[27] 


185 


XTE J0929-314 


[12] 


45 


EXO 0748-676 


[3^ 




FIGURE 1. The K-S probability distribution of the synthetic populations used to sample the joint 
P — P parameter space. The relaxed 2D K-S filter allows additional acceptance of populations with 
0.05 < P2DK-S < 0.2 by oversampling the extreme outliers of the P P distribution in order to probe 
for possible contaminations and extreme fluctuations. The distribution also shows how optimally the 
P — P parameter space is sampled with a peak sampling rate around the nominal acceptance value of 
PiDK-s ~ 0.2. The dotted line is the relaxed 2D K-S filter for the synthetic populations that is used to 
construct the predictive distribution in Fig 2. The solid line is the conventional K-S filter that would only 
accept strictly consistent P P samples.^ 



^ Labelled version of Fig. 1 in ll2Tll 




(a) The P P distribution of millisecond radio pulsars evolved 
from a progenitor population consistent with the observed mil- 
Usecond X-ray pulsars (Table 1). The probability map is a pro- 
duced by seeds sampled from a discrete progenitor millisecond 
X-ray pulsar population.^ 



(b) In this case, the probability map is a produced by seeds sam- 
pled from a continuous millisecond X-ray pulsar progenitor pop- 
ulation consistent with the observed milUsecond X-ray pulsars 
(Table 1). 



FIGURE 2. The P — P distribution of millisecond pulsars for the standard model vs. observed mil- 
lisecond radio pulsars. The color contours are the expected millisecond radio pulsar distribution for the 
observed millisecond X-ray pulsars, with red representing the highest density. Any millisecond radio 
pulsar outside of the color contours cannot be produced by the observed millisecond X-ray pulsars with 
more than 95 % confidence. The dashed red line is the 68 % confidence limit of the expected P — P 
distribution taking the observed millisecond X-ray pulsars as the progenitor population. Millisecond radio 
pulsars within the shaded region may be produced by millisecond X-ray pulsars with spin distributions dif- 
ferent than what is observed. The area shaded with lines assumes a maximum spin frequency Vmax = 619 
Hz, which the fastest spinning observed millisecond X-ray pulsar (i.e. 4U 1608—52). The shaded area 
extends to the dotted region if the maximum spin frequency for millisecond X-ray pulsars is allowed to be 
Vmax — 760 Hz which is the theoretical upper limit predicted by [8J. The blue squares and red filled circles 
are the observed milUsecond radio pulsars in single and binary systems respectively. The area outside of 
the shaded region is not sensitive to the prior, i.e. the observed millisecond radio pulsars outside of the 
shaded area cannot be produced consistently by the standard model for any millisecond X-ray pulsar spin 
distribution with more than 95 % confidence. The spin-down values for the observed milUsecond radio 
pulsars are corrected for secular acceleration ll33l l6ll. The spin-up line (P ^ iiiPq^^) for m = MecUI and the 
characteristic age line for = lO'" yrs are shown with dashed and dash-dotted lines. ^ 



^ Fig. 2 in ll2lll optimized for grayscale. 

^ High resolution full color figures available at URL: I http : // w w w.kiziltan . org/research.htrnI| 



